DESY 99-099 
MPI/PhT/99-029 



|hep-ph/9907489 



July 1999 



ISSN 0418-9833 



Ultrasoft Effects in Heavy-Quarkonium Physics 

Bernd a. Kniehl* 

Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut), 
Fohringer Ring 6, 80805 Munich, Germany 

Alexander A. Penin^ 

II. Institut fiir Theoretische Physik, Universitat Hamburg, 
Luruper Chaussee 149, 22761 Hamburg, Germany 



Abstract 

In the fi'amework of nonrelativistic QCD, we consider a new class of radiative cor- 
rections, which are generated at next-to-next-to-next-to-leading order through the 
chromoelectric dipole interaction of heavy quarkonium with ultrasoft virtual gluons. 
We provide analytical formulae from which the resulting shifts in the quarkonium 
energy levels and the wave functions at the origin may be calculated. We discuss 
the phenomenological implications for the top-antitop and T systems and point out 
some limitations of describing charmonium using a Coulomb potential. 
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1 Introduction 



Recently, essential progress has been made in the theoretical investigation of the pair 
production of heavy quarks at threshold, and analytical results are now available up 
to next-to-next-to-leading order (NNLO) |l]J|J|,|,|,|,0g . The NNLO corrections have 
turned out to be so sizeable that it appears to be indispensible to also gain control over 
the next-to-next-to-next-to- leading order (N^LO), both in regard of phenomenological 
applications and in order to understand the structure and the peculiarities of the threshold 
expansion. In this paper, we take the first step in this direction and investigate a particular 
class of N^LO corrections, namely those arising from the emission and absorbtion of virtual 
ultrasoft gluons by the heavy quarks. Such ultrasoft corrections are absent in lower orders 
and represent a genuinely new feature of the N^LO. 

Specifically, we study the near-threshold behavior of the vacuum-polarization function 
n(g^) of a heavy-quark vector current = qjfj.q, 

{q,qu - 9,uq') = ^ J d'xe^'^ ^0\Tj,{x)jM\0) ■ (1) 

Its imaginary part gives (up to a constant factor) the spectral density of qq production in 
e~^e~ annihilation. Near threshold, the heavy quarks are nonrelativistic, so that one may 
consider the quark velocity /3 (or inverse quark mass) as a small parameter. An expansion 
in (3 may be performed directly in the Lagrangian of quantum chromo dynamics (QCD) 
by using the framework of effective field theory. In the threshold problem, there are four 
different scales (i) the hard scale (energy and momentum scale like m^); (ii) the 
soft scale (energy and momentum scale like /3mq); (iii) the potential scale (energy scales 
like P'^mg, while momentum scales like I3mq); and (iv) the ultrasoft scale (energy and 
momentum scale like (3'^mq). The ultrasoft scale is only relevant for gluons. 

By integrating out the hard scale of QCD, one arrives at the effective theory of non- 
relativistic QCD (NRQCD) |T0[. If one also integrates out the soft scale and the potential 
gluons, one obtains the effective theory of potential NRQCD (pNRQCD), which contains 
potential quarks and ultrasoft gluons as active particles . The dynamics of the quarks 
is governed by the effective, nonrelativistic Schrodinger equation and by their interac- 
tion with the ultrasoft gluons. To get a regular perturbative expansion within pNRQCD, 
this interaction should be expanded in multipoles. The corrections from harder scales are 
contained in the Wilson coefficients, leading to an expansion in «<;, as well as in the higher- 
dimensional operators of the nonrelativistic Hamiltonian, corresponding to an expansion 

in l/lTLq. 

This paper is organized as follows. In Section 2, we introduce the basic features of 
the pNRQCD formalism and derive formulas from which the leading ultrasoft corrections 
to the energy levels and the wave functions at the origin of heavy quarkonia may be 
evaluated. In Section 3, we present a numerical analysis and discuss phenomenological 
implications of our results. Section 4 contains our conclusions. In the Appendix, we 
explain how the "QCD Bethe logarithms" introduced in Section 2 may be reduced to 
one-dimensional integrals of elementary functions. 
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2 Ultrasoft corrections to energy levels and wave func- 
tions 



In this section, we briefly recall the formalism of pNRQCD and calculate the leading 
ultrasoft corrections to the energy levels and the wave functions at the origin of heavy 
quarkonia. The basic quantity of pNRQCD is the nonrelativistic Green function, G*'", of 
the Schrodinger equation, 

{W^^ - E) G^-°(x, y , E) = 5(3) (x - y) , (2) 

where W'" is the nonrelativistic Hamiltonian for the quark pair in the colour- singlet (s) 
[colour-octet (o)] state, defined by 

W^" = -— + V''\x) + . . . , 
nig 

V^^'°(x) = V^^\x) + . . . , (3) 

with Ax = and x = |x|. The ellipses stand for the higher-order terms in as and 1/mq. 
The Coulomb (C) potentials for the singlet and octet states are attractive and repulsive, 
respectively, and are given by 

VSix) = -Cf^, 

X 

Vi(x) = (§i - C) ^, (4) 

where Ca = 3 and Cp = 4/3 are the eigenvalues of the quadratic Casimir operators of 
the adjoint and fundamental representations of the colour group, respectively. 

The propagation of the quark-antiquark pair in the singlet and octet states is described 
by the singlet and octet Green functions, respectively, which have the following spectral 



representations 12 



^ ^ ^1 En-E Jo 



(27r)3 kynig-E ' 

^^""'^'^^-io (2vr)3 kymg - E ' 

where ipm and ip^" are the wave functions of the qq bound and continuum states. Note 
that a discrete part of the spectrum (bound states) only exists for the singlet Green 
function. 



The nonrelativistic expansion in as and (3 = y 1 — Arriq/s, where ^/s is the qq centre- 
of-mass energy, provides us with the following representation of the heavy-quark vacuum- 
polarization function near threshold: 

U{E) = ^Gh{0,0,E) + ..., (6) 
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where E = ^/s — 2mq is the qq energy counted from the threshold and is the leading- 
order Coulomb Green function, which sums up the (a^/ /?)" terms singular near the thresh- 
old. The ellipsis stands for the higher-order terms in as and (3. 

We are interested in the correction to the Green function induced by a virtual ultrasoft 
gluon. The corresponding diagram is shown in Fig. 1. The leading contribution is the 
one due to the chromoelectric dipole interaction ^'^(rq — rg) • E of the quark- antiquark 
pair with the ultrasoft gluon |jl2|. We wish to calculate the corresponding corrections 
to the energy levels and the wave functions at the origin of several low-lying resonances, 
which represent key parameters for the analysis of the threshold production of top- and 
bottom-quark pairs. Near the nth pole of the discrete spectrum, the correction to Gq 
reads 



where 



J{E) = C,gl ^(r.)k.(r,)„ur^(i? - /m,) , (8) 

with 

in D space-time dimensions. After integration over /q? we recover the well-known non- 
relativistic perturbation theory expression. The pole in the nonrelativistic propagator is 
bypassed according to the standard prescription E —>■ E+ie. The remaining integral over 1 
is ultraviolet divergent for = 4. To obtain a finite result, we use dimensional regulariza- 
tion with D = 4 — e. The nonrelativistic perturbation theory of quantum electrodynamics 
(QED) in dimensional regularization is comprehensively described in Refs. JTSIJT^ . We 
thus obtain 

/^^■(p)=p'^f- + ln^ + ^-ln2') , (10) 

where 1/e = 1/e + [ln(47r) — 7£;]/2. Note that this divergence is spurious. It arises in 
the process of scale separation due to the use of pNRQCD perturbation theory at short 
distances where it is inapplicable. In the total N^LO result, the pole in 1/e is cancelled by 
the infrared poles coming from the hard- and soft-scale corrections. However, since these 
corrections are still unknown, we subtract the divergent part according to the MS scheme. 
This means that the same scheme must be used for the calculation of the hard- and 
soft-scale corrections. As a consequence, the partial result for the ultrasoft contribution 
depends on the auxiliary "factorization scale" fif, which drops out in the total result. 

The matrix element (r)k„ is the one between the singlet Coulomb wave function of 
principal quantum number n and the octet Coulomb wave function of momentum k. Since 
this matrix element is finite, we may use the wave functions of the Schrodinger equation 
in three dimensions, with e = 0. Writing 

En = + ^En , 
\MO)\' = \^^iO)Wl + Arn) , (11) 
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where the Coulomb values are 



E. 



c 



with Xs = asCFmq/2, we obtain the leading ultrasoft corrections as 



dJ{E) 



dE 



E=EC 



121 



(13) 



Inserting Eq. (H), we thus obtain 



2CFas 
3^ Jq (27r)3 



\{v)^n?{E^;-}^'/m^ 



In 



En ~ k'^/rnc^ 



(14) 



Except for the terms involving ln[£'p/(£'^ — /c^/mg)], we may evaluate the right-hand 
sides of Eq. ([T^) analytically. In fact, making use of the completeness relation, 



d^k 
(27r)3 



(15) 



this problem may be reduced to the calculation of the diagonal matrix elements of appro- 
priate local operators. 



rlE^-nn''r 



{e^ - 



Vy- _ v°f + —(V'- V") (e^ - 3V' + V) - 2?^" 

m.. V / 



c 



niq 

-C^ + -C^Cf + ; 

4 n \n 



m 



1 I nn 



n 



r\V' - V°f + — (e^- 2V' + V° 



1 



(16) 



The contribution is purely Abelian and coincides with the QED result for the positron- 
ium bound state |jl3| . The "maximal non- Abelian" contribution proportional to C\ in the 
local part of the energy shift may be read off from the analysis of the infrared properties of 
the non- Abelian static potential UTS). The C^/e term in the ultrasoft contribution to the 
energy shift cancels the corresponding infrared pole in the potential |l^. Note that there 
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is no such term in the wave-function correction. This may be understood by observing 
that, in contrast to the case of the energy shift, the wave function receives an additional 
infrared C\/ e contribution from the hard matching coefficient of the heavy-quark vector 
current ||T^. This hard contribution must cancel the potential-related divergence in the 
correction to the wave function, but it does not affect the energy levels. 

The logarithmic contributions in Eq. (|l^) represent a pure "retardation effect" and 
cannot be reduced to local-operator contributions. It is convenient to introduce the "QCD 
Bethe logarithms". 



1 



CpE^ 



rV' _ 



(27r)3 
\{rhn\' 



e: 



c 



k'^/nig] In 



c 



E^ - k'^/ruq 



(17) 



In the Appendix, we explain how these QCD Bethe logarithms may be reduced to one- 
dimensional integrals of elementary functions. For n = 1,2,3, we obtain the following 
numerical values: 



Lf = -81.5379^ 
Lf = -5.7675, 



Lf 



-37.6710, 
0.7340, 



tE 
^3 

-^3 



-22.4818, 
2.2326. 



The final results for the ultrasoft corrections to the heavy-quarkonium energy levels and 
wave functions at the origin read 




1 



TT 



E? 
1 
4 



C\Cf H — 7;CaC% H — -C^ 



n 




\n2\+C%Lt 



(19) 



3 Phenomenological applications 

We are now in a position to present a numerical analysis and to discuss the phenomeno- 
logical implications of our results. As illustrative examples, we consider, in Figs. 2 and 
3, the bb and tt ground states, with n = 1, respectively, and investigate the size and the 
dependence on the factorization scale /^/ of the ultrasoft corrections to the energy level 
El and the square of the wave function at the origin |'?/'i(0)p. As input values for the 
pole masses and the strong coupling constant, we use mi, = 4.8 GeV, rrit = 175 GeV 
and as{Mz) = 0.118, respectively. Note that, in Eq. (p!9|), one power of as refers to the 
ultrasoft gluon interaction and should be evaluated at the ultrasoft scale, a^m^, while 
the two residual powers of as originate from the Coulomb Green function and should be 
evaluated at the Coulomb scale, asTUg. This leads us to solve the functional equations 
asia'^nig) = as and as{asmq) = as for the ultrasoft and Coulomb regimes, respectively. 
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In the top-quark case, we have as = 0.146 at the Coulomb scale and as = 0.196 at the 
ultrasoft scale. Both solutions are in the perturbative region. In the bottom-quark case, 
we obtain as = 0.34 at the Coulomb scale and as = 0.47 at the ultrasoft scale. The last 
value seems to be too large for a reliable perturbative calculation. Thus, to be on the safe 
side, we redefine the ultrasoft scale to be Sa'^rub. This leads to as = 0.34, which coincides 
with the value at the Coulomb scale. The fif dependence is cancelled in the total result 
by hard terms of the form ln{fif/mq) and by soft terms of the form ln{fif /asiriq). In order 
to minimize these unknown logarithmic contributions and to absorb the large logarithms 
into the know ultrasoft terms, we select Hf from the interval asfriq < fif < niq. 

From Figs. 2a and 3a, we observe that, depending on /ij, the ultrasoft N^LO corrections 
to the energy levels of the bb and ti systems may be as large as —120 MeV and +120 MeV, 
respectively. As we see from Figs. 2b and 3b, the ultrasoft N^LO wave-function correction 
reaches —40% for the bb system and —7% for the tt system. The circumstance that the 
corrections are close to zero at some points of the interval a^mg < fif < rrig justifies this 
choice of factorization scale. It is interesting to compare the ultrasoft N^LO corrections 
presented in Figs. 2 and 3 with the corresponding NLO and NNLO corrections, which we 
extract from Ref. 0. In the bottom-quark case, the NLO (NNLO) contributions to AEi 
and Aipf approximately amount to -100 MeV (-200 MeV) and +50% (+150%), respec- 
tively. In the top-quark case, the corresponding values are —700 MeV (—700 MeV) and 
—15% (+15%), respectively. This comparison nicely illustrates the numerical significance 
of the new N'^LO corrections. 

Interesting, relevant and timely applications of our results include the studies of top- 
quark pair production at threshold, of low-lying T resonances and of T sum rules. In the 
case of the IS T resonance, where local duality is expected to work, it is interesting to 
compare the perturbative ultrasoft contribution, related to the scale a^m^ ^ 1 GeV, with 
the leading nonperturbative contribution of the gluon condensate, due to nonperturbative 
fluctuations at the scale Aqcd, which originates from a diagram of the type shown in Fig. 1 
with a broken gluon propagator. In the case of the IS energy level, the latter is given by 



Using (asG^) = 0.06 GeV, we have AEi ^ 60 MeV, which is comparable to the ultrasoft 
contribution. 

In the case of the J/ if) resonance, which is sometimes optimistically considered to be 
a Coulomb system, the ultrasoft scale is of the order of Aqcd if one assumes ag to freeze 
at 1 GeV. Thus, one has to accurately separate perturbative and nonperturbative con- 
tributions. Nevertheless, our results suggest that the J/ip resonance cannot be regarded 
as a Coulomb potential system. Indeed, the potential model can be destroyed either by 
the large contribution from the gluon condensate, which is nonpotential because the non- 
perturbative scale Aqcd is below the potential scale a^mc ~ 700 MeV, or by the large 
perturbative retardation contribution. The first one is proportional to l/ctf, while the 
second one is proportional to a^. We observe that, due to their different dependence on 
as, these two contributions cannot be simultaneously small for any low-scale evolution 




(20) 
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of the coupling constant. For rric ~ 1.4 GeV and as ~ 0.5, the typical size of nonpo- 
tential contributions to the J/ip energy levels is about 400 MeV, which is comparable 
to the inverse Bohr radius. This clearly indicates that the J/ip system is far from being 
Coulombic. 



4 Conclusions 

We conclude with a few general remarks concerning the structure of higher-order cor- 
rections in pNRQCD. At NLO, the only source of corrections is the running of as- At 
NNLO, higher-dimensional operators start to contribute. At N^LO, retardation effects, 
which cannot be described by local operators, enter the game. The leading retardation ef- 
fects, which are under consideration here, arise from the chromoelectric dipole interaction 
of heavy quarkonium with virtual ultrasoft gluons as depicted in Fig. 1. To our knowledge, 
these effects have not been studied elsewhere in the literature. We emphasize that they 
constitute a genuinely new feature, which is absent in NLO and NNLO. In particular, 
they are not contained in any of the popular renormalon-based mass definitions. On the 
other hand, they are expected to be the last source of unexpectedly large corrections. 
Thus, our result sets the scale of the N^LO corrections. 

In this paper, we took a first step towards the N^LO analysis of the heavy-quark 
threshold dynamics. Our results entirely account for the ultrasoft-scale physics and should 
be complemented by the soft and hard contributions. The analysis of the last one includes 
the calculations of the three-loop hard matching coefficient and the three-loop potential. 
The corresponding two-loop results may be found in Refs. |ll7| , p!8| , respectively. Some 



N^LO results on the effective Lagrangian were presented in Ref. [ITU 
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Appendix 



Only the S-wave component of the (colour-singlet) Green function, with angular-momen- 
tum quantum number / = 0, contributes to its value at the origin and, therefore, to 
the leading ultrasoft corrections to the vacuum-polarization function. This means that 
only the I = 1 component of the colour-octet wave function contributes to the matrix 
element (r)kn. The corresponding Coulomb wave functions for the attractive (singlet) 
and repulsive (octet) potentials read 

i'noir) = -^=Rno{r) , 
v47r 



«,(r)=e-<'/-^f)^fl„(,)!i-i:, (21) 



n 



where 6^ is the / = 1 Coulomb phase and 

fa— / 2,1 \ 1/2 

Rkiir) = ^kr f e*''^F(2 + zu, 4, -z2kr) . (22) 

Here, F is the confluent hypergeometric function, u = X^pi/k and 

I f Ca \ I 
Pn = -h^-l =-. 23 
n \2Cf / on 

After some algebra, we obtain the QCD Bethe logarithms of Eq. ( p!?]) in terms of one- 
parameter integrals of elementary functions, as 

/•oo 

L^'^ = / duY^^'Hu)X'^{u) , (24) 
Jo 



where 



and 



Xi(z/)=pi + 2, 

I/2(2p2 + 9p2 + 8)-p2(p2+4) 



X2M 



u\Spl + 60p2 + 123p3 + 66) - 2z/2p2(6p2 + Alp, + 54) + 3p^(p3 + 6) 

w^+w ■ ^ ^ 

Since it is usually sufficient to consider n = 1,2,3 in practical applications and the 
expressions for X„ with n > 3 are somewhat cumbersome, we refrain to listing the latter. 



11 



Figure captions 



Fig. 1. Feynman diagram giving rise to the ultrasoft contribution at N^LO. The single 
and double lines stand for the singlet and octet Green functions, respectively, the wavy 
line represents the ultrasoft-gluon propagator in the Coulomb gauge, and the vertices 
correspond to the chromoelectric dipole interaction. 

Fig. 2. Ultrasoft corrections to (a) the energy level Ei and (b) the square of the wave 
function at the origin |'?/'i(0)p as functions of the factorization scale fij for the bb ground 
state, with principal quantum number n = 1. 

Fig. 3. Same as in Fig. 2, but for ti. 




Fig. 1 



12 




/X/ (GeV) 

Fig. 2a 



Aijj X 10^ -25 




fif (GeV) 
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Fig. 3b 
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